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A  STUDY  OF  PLASMA  DECAY  IN 
A  TOROIDAL  MAGNETIC  FIELD 
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and  A.  P.  Zhilinskiy 

Leningrad  Polytechnical 
Institute  iraeni  M.  I.  Kalinin 


This  article  presents  the  results  of  a  theoretical 
and  experimental  study  of  the  decay  of  slightly 
ionized  plasma  in  a  toroidal  chamber  along  which  a 
magnetic  field  is  directed.  Equations  for  plasma 
decay  have  been  derived  from  the  theoretical 
study  and  their  solution  found  and  analyzed  as  a 
description  of  the  spatial  distribution  of  the 
plasma  concentration  and  of  the  changes  in  its 
concentration  over  time.  The  change  over  time  in 
the  concentration  of  charged  plasma  particles 
undergoing  decay  within  the  helium-filled  dielec¬ 
tric  chamber  was  calculated  experimentally  for  a 
variety  of  conditions.  Experimental  data  on  the 
effect  of  toroidal  drift  in  the  magnetic  field  on 
plasma  decay  agree  with  the  theoretical  conclu¬ 
sions. 


Introduction 

M#st  of  the  research  in  the  area  of  plasma  decay  within  a  magnetic 
field  is  directed  toward  the  end  cf  obtaining  data  on  charged  par¬ 
ticle  diffusion.  In  one  of  their  first  experiments,  Bostick  and 
Levine  [1]  studied  the  decay  of  slightly  ionized  plasma  within  a 
metal  toroidal  chamber.  This  work  contains  their  calculation  of 
the  decay  time  constant  based  on  curves  describing  the  change  over 
time  in  the  concentration  of  the  charged  plasma  particles.  The 
relationship  of  the  decay  constant,  as  it  turns  out,  to  the  longi¬ 
tudinal  magnetic  field  has  a  maximum,  the  value  of  the  maximum 
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decay  constant  increasing  with  the  pressure  of  the  gas.  In  the 
authors'  view,  the  presence  of  this  maximum  provides  evidence  of 
an  irregularity  in  the  relationship  of  the  coefficient  of  diffusion 
to  the  magnetic  field  as  predicted  by  the  theory. 

Work  [1],  as  well  as  many  subsequently  published  articles,  considers 
the  experimental  results  referred  to  as  proof  of  the  occurrence  of 
plasma  instability  as  well  as  of  an  anomalous  diffusion  accompanying 
it.  It  should  be  pointed  out,  however,  that  this  is  a  questionable 
interpretation  of  the  experimental  data,  considering  the  fact  that 
there  had  been  no  analysis  of  the  effect  on  the  plasma  decay  of  the 
charged  particle  drift  caused  by  the  nonuniformity  of  the  magnetic 
field  (although  the  authors  did  refer  to  the  possibility  of  a  sub¬ 
stantial  effect  of  toroidal  drift).  The  purpose  of  the  present 
work  has  been  precisely  to  study  the  effect  on  plasma  decay  of 
charged  particle  drift  in  a  nonuniform  magnetic  field. 

1.  THE  THEORY  OF  PLASMA  DECAY  IN  A  TOROIDAL  MAGNETIC  FIELD 

a)  The  Decay  Equation 

Let  us  consider  the  decay  of  slightly  ionized  plasma  in  a  toroidal 
chamber  along  which  a  magnetic  field  is  directed  (Figure  1).  With¬ 
in  the  cylindrical  coordinate  system, 
the  field  components  are  calculated  in 
accordance  with  the  eolations 

«;=//, //,=//.= o.  (1) 

In  the  course  of  our  discussion  we 
will  consider  the  radius  of  the  cham¬ 
ber's  curvature  to  be  much  greater 
than  its  transverse  dimensions 


Figure  1  *>L.  (2) 


s 


We  calculate  the  change  in  concentration  of  charged  particles 
during  the  process  of  plasma  decay  in  accordance  with  the  equation 


at 


■^r(n.».)  =  0, 


(5) 


in  which  n.  -  concentration,  «.  -  average  particle  velocity;*  we 
assume  that  the  processas  of  the  removal  of  charged  particles  en 
masse  are  immaterial.  The  condition  of  quasineutrality  of  a  plasma 


*Index  *  =  «  will  indicate  electron  values;  index  «=-/—  the  ion  values. 
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consisting  of  electrons  and  singly-charged  ions  results  in  the 
equality  of  their  concentrations 

n,  =  H{  =  n, 

(*> 


and,  correspondingly,  to  equality  in  the  divergence  of  the  flows 

V  (*,«.)  =  V(n<«,).  ^ 

Average  charged  particle  velocity  may  be  calculated  employing  the 
averaged  equation  for  particle  motion,  which,  as  we  know,  is  a 
corollary  of  the  kinetic  equation 


+ 


m. 


it 

X  H] 


—  Z.eE- 

T.9n . 


-m.v  _ 


(6) 


Here  m.  -  mass;  Z.e  -  charge  (Z«  =  1,  Z.  =  — 1)  ;  e  -  electrical  field 

intensity;  H  -  magnetic  field  intensity,  7*.  -  temperature  in  energy 
units  (temperature  assumed  to  remain  constant);  ».  -  effective  fre¬ 
quency  of  the  collisions  of  *-type  particles  with  neutral  atoms 
(the  frequency  of  the  collisions  of  charged  particles  with  neutral 
atoms  in  a  slightly  ionized  gas  is  much  greater  than  the  frequency 
of  the  collisions  of  charged  particles  with  one  another). 


Average  particle  velocity  is  usually  much  lower  than  thermal  veloc¬ 
ity  ^  ^  |hT7  ,  and  we  may  ignore  the  second,  nonlinear  term  in 

equation  (6).  In  the  case  of  a  sufficiently  slow  plasma  decay, 

in  which  >  we  may  ignore  the  first  term  as  well.  Equation 

(6)  takes  the  form 


- ^[«.XHJ  =  Z.#E  -  T.^-.  (7) 

By  solving  the  vector  equation  we  obtain  the  following  expressions 
for  the  transverse  (perpendicular  to  the  magnetic  field)  components 
of  the  average  velocity  * 


+  The  longitudinal  component  of  average  particle  velocity  «.f  in  the 
case  of  a  closed  torus  is  equal  to  zero. 
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_MEXHi 

■ia  /  *i 


<  H]  cTa  [Vn.  X  H] 


_  Z««*,E» _ '<aT,Vln, 


—  itL 

in  which  ™,c  -  Larmor  frequency;  index  t  indicates  the 

transverse  vector  components. 

The  vector  represents  the  velocity  of  particle  drift  conserved 
in  the  absence  of  collisions.  Let  us  note  that  the  drift  produced 
by  the  nonuniform  magnetic  field  is  implicitly  given  by  the  second 
in  the  expression  for  uu  [see  (9)1.  Vector  ua  represents  the  ve¬ 
locity  of  particle  movement  along  the  electrical  field  and  the 
concentration  gradient  connected  with  the  collisions. 

Taking  into  consideration  (1)  and  (2),  we  may  easily  obtain  by  em¬ 
ploying  (8)  the  following  expression  for  the  divergence  in  particle 
flow 


* («.*.) = ~  \-rn-^rr'  - 

%  _L  J_  f  1 7V",  1  d  r  r  ‘  dn,  1 

*  r  dr  [  (*>!  ^  J  m*  r  0r  [  (*’  -h  -1)  °r  | 


Za<-v.  1  d 


o  (n,Er)  Z.*v. 

0*  '".(^  +  ''2) 

dz 

*«.  .  2 cT. 

!)  /  vJ  ' 

j*  d* 

The  last  term  defines  the  change  in  concentration  due  to  the  drift 
caused  by  the  nonuniform  magnetic  field. 

In  order  to  state  expression  (8)  for  flow  divergence  more  precisely, 
we  must  calculate  the  intensity  of  the  electrical  field  in  the 
plasma  from  quasineutrality  conditions  (4)  and  (5).  We  will  keep 


in  mind  in  our  calculations  that  the  following  inequalities  are 
satisfied  for  the  experimental  conditions  discussed  below 


<B*  ^  W«I  ">,<*>(  ^ 

(the  first  inequality  being  a  corollary  of  the  second). 


(10) 


Let  us  observe  that  in  the  case  of  plasma  decay  within  a  straight 
cylinder  of  great  length  (when  R  -*co),  the  quasineutrality  condi¬ 
tion  (5)  may  be  satisfied  if  the  components  of  electron  and  ion 
velocity  in  the  direction  of  the  gradient  of  concentration  are  the 
same  (u,x=iijA)  .  In  this  instance  we  obtain  a  regime  of  ambipolar 
diffusion.  While  maintaining  inequality  (10),  the  rate  of  diffusion 
yields  the  relationship  (see,  for  example,  [2]) 


n  —  (y.-t-  7\)v, 

•L  ~  2 


(ii) 


The  electrical  field  insuring  the  quasineutrality  of  the  plasma 
is  calculated  employing  the  following  inequality 


T.V,n 

en 


(12) 


The  electrical  field  within  the  toroidal  chamber  differs,  of  course, 
from  (12).  We  will  look  for  it  in  the  form 


E  =  E2>-+-E8.  (13) 

Employing  (12)  and  (13)  and  equating  the  divergences  in  the  flow 
of  electrons  and  ions  (5),  (9)*  we  obtain  the  following  equation 
after  simple  transformations 


1  t 

r  Or 


:nE„ 


re*  .nEt 


//(o 


I) 

t't.nEu 


».(-?- 


vJ)K  i  {' 


cv*n£. 


«.(-!■ 


U) 


2cn  {T ,  T j)  \  = 
eRH 


j=° 


(W) 


[inequalities  (2)  and  (10)  here  taken  into  consideration]. 

We  may  calculate  the  components  of  the  electrical  field  separately 
by  setting  the  brackets  under  the  derivative  sign  equal  to  zero 


f  —1±Ll±J*L.  r  _  . 

£*'  ^  —  ^  ’ 

R  =  E  -+-E  — Z*.lUL _ 3  (7*,  Tj)  yH  2 (T,  Tt)  [H  x  yH] 

°  *  *  n  •  ir  «v(  ift 


(15) 
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It  should  be  observed  that  the  distribution  of  the  space  charge 
field  calculated  in  this  manner  satisfies  not  only  equation  (14), 
but  also  the  condition  of  the  equality  of  electron  and  ion  flows 
at  any  point  along  the  boundary  of  the  plasma  (at  0  boundary  con¬ 
centration).  This  condition  is  obviously  necessary  in  the  case 
of  dielectric  chamber  walls.  It  is  possible  that  the  field  dis¬ 
tribution  in  (15)  is  not  the  only  one  satisfying  equation  (14). 

Substituting  (15)  in  expression  (9)  for  V(n,».)  and  keeping  (2)  and 
(10)  in  mind,  we  obtain  the  following  relationship 


V  (im.)  =  V  (nu<)  =  —  ,  (16) 

in  which  are  introduced  the  designations 


d*n 
dt*  ' 


tv 


2^, 

R 


— rate  of  toroidal  plasma  drift. 


n  --  Tj) 


(17) 

(18) 

(19) 


— coefficient  of  ambipolar  diffusion  in  the  absence  of  the  magnetic 
field  (or  along  the  field). 


The  equation  for  plasma  decay  (3)  thus  takes  the  form 


dn 

Hi 


— D 


dn  n 


(20) 


b)  Integrating  the  Decay  Equation 


Assuming  concentration  on  the  plasma  boundary  (on  the  surface  of 
the  toroidal  chamber)  is  equal  to  zero,  let  us  now  find  the  solu¬ 
tion  to  the  equation  for  plasma  decay  (20).  The  general  solution 
to  the  equation  may  be  presented  in  the  form 


a 


—  2***  'kst(r,  z)eiDl 


I 

(21) 


where 


J _ .  «UJ 

**  A*  40x  .  (22) 
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■**  -  eigenfunctions  of  the  equation 


satisfying  the  zero  boundary  conditions;  At  -  the  corresponding 
eigenvalues;  ck  -  constant  coefficients  designating  the  initial 
distribution  of  concentration. 

After  the  passage  of  a  sufficient  period  of  time  following  the 

beginning  of  plasma  decay  *^257  only  the  first  term  in  sum  (21), 

corresponding  to  the  greatest  value  a  =  A0,  proves  to  be  of  sub¬ 
stance 


-  v  .  .  u>. r 

n  =  c o*  s0(r,  z)e  x  . 


Plasma  decay  in  this  stage  proves  to  be  exponential.  In  accordance 
with  (22),  the  inverse  constant  of  decay  may  be  written  in  the  form 
of  the  sum  of  the  term  defining  the  rate  of  diffusion  elimination 
of  the  particles  (unconnected  with  the  curvature  of  the  chamber) 
and  the  term  appearing  as  a  result  of  toroidal  drift 


_L  =  .  _L  Jii  -  _L  ( 

A  l  4  Dx  -  2t„  \D0  ■+■  )  ’ 

JL  —  JL. _ 2D II  n  —  mA0 _ °IIA0 

A°  RA'  ’  #  2  *'  ’  (25) 

As  can  be  seen  from  (25),  with  a  monotonic  change  in  Dj.  ,  the  time 
constant  passes  through  its  maximum  when  the  coefficient  of 
transverse  diffusion  equal  to  D0  , 

The  function  s0(r,  z)  entering  into  (24)  defines  in  accordance  with 
(25)  what  is  referred  to  as  the  "diffusion”  distribution  of  con¬ 
centrations  over  a  cross  section  of  the  chamber.  We  obtain  this 
distribution  with  the  "rectification"  of  the  toroidal  chamber 
(when  R  -*  qo  ),  in  the  case  of  a  toroidal  chamber  with  a  finite 
radius  of  curvature,  the  distribution  maximum  appears  in  a  position 
shifted. farther  in  the  direction  of  the  drift  [the  exponential 
factor  in  (24)]  the  greater  is  the  ratio  of  the  rate  of  drift  to 
the  coefficient  of  lateral  diffusion. 

The  concentration  distribution  in  the  initial  stage  of  plasma  de¬ 
cay  may  differ  from  (24).  In  this  instance,  the  rate  of  the  elim¬ 
ination  of  charged  plasma  particles  may  be  either  higher  or  lower 
than  at  the  end  of  the  decay  process. 


iii*2 


vie*' 


We  should  point  out  that  the  drift  of  the  charged  particles  to¬ 
ward  the  chamber  walls  occurs  as  a  result  of  lateral  diffusion 
alone  and  is  determined  entirely  by  the  concentration  gradient. 

The  particle  flow  connected  with  the  drift  nv>  in  fact  vanishes 
near  the  walls,  since  the  concentration  vanishes  as  well.  If  the 
initial  distribution  of  concentrations  is  symmetrical  (near  dif¬ 
fusion  [concentration — TR]),  the  rate  of  particle  drift  toward  the 
walls  therefore  proves  to  be  identical  to  that  which  occurs  in  the 
absence  of  toroidal  drift.  In  consequence  of  the  drift  there  then 
occur  a  shift  in  the  distribution  maximum  and  increases  in  both 
the  rate  of  plasma  diffusion  and  the  concentration  gradient  near 
the  chamber  walls,  toward  which  the  drift  is  occurring.  Practi¬ 
cally  exponential  decay  with  the  time  constant  calculated  in  ac¬ 
cordance  with  equality  (25)  must  be  established  over  a  period  of 
time  of  the  order  of  the  "drift  time"  /Aj  . 

c)  Comparing  Theoretical  with  Experimental  Results 

Comparison  of  theoretical  conclusions  with  the  experimental  re¬ 
sults  obtained  by  Bostick  and  Levine  [1]  shows  a  good  correspondence 
between  them.  The  maximum  for  the  curve  defining  the  relationship 
between  the  plasma  decay  constant  and  the  magnetic  field  -(H)  ob¬ 
served  in  work  [1]  probably  corresponds  to  the  maximum  following 
from  relationship  (25). 

The  maximum  value  of  the  decay  constant  ,  calculated  on  the 
basis  of  the  experimental  data,  increases  with  pressure,  which  also 
agrees  with  the  formula  in  (25).  We  should  point  out  that  there  is 
a  substantial  quantitative  discrepancy  in  value  .  The  value 
determined  experimentally  proves  to  be  1.5-2. 5  times  smaller  than 
the  value  as  calculated. 

Quantitative  agreement  between  the  conclusions  drawn  from  the 
theory  developed  above  and  the  experimental  data  of  Bostick  and 
Levine  was  hardly  to  be  anticipated.  These  experiments  were  con¬ 
ducted  within  a  metal  chamber.  Equalizing  the  potential  on  the 
boundaries  of  the  plasma  should  lead  to  a  change  in  the  space- 
charge  field  calculated  in  accordance  with  relationships  (15)  for 
any  case  near  the  boundaries.  The  comparatively  large  curvature 
of  the  toroidal  chamber  employed  in  the  experiments  could  be  an¬ 
other  cause  of  discrepancies.  The  toroid's  radius  of  curvature 
was  altogether  2.5  times  greater  than  the  width  of  the  rectangular 
cross  section,  and  the  low  toroidality  condition  (2)  was  not  sat¬ 
isfied.  Discrepancies  between  the  experiment  and  the  theory  may 
also  have  been  connected  with  deficiencies  in  the  procedure  used 
in  work  [1]  to  measure  concentration.  The  superhigh-frequency 
measuring  field  employed  in  the  experiments  could  cause  string 
plasma  disturbances,  since  the  measurments  were  conducted  at 
frequencies  near  that  of  electronic  cyclotron  resonance. 

In  order  to  test  the  theory  quantitatively,  we  undertook  an  experi¬ 
mental  study  of  plasma  decay  within  a  dielectric  toroidal  chamber 
with  a  radius  of  curvature  much  greater  than  its  transverse  dimen¬ 
sions. 
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2.  THE  METHOD  OF  THE  EXPERIMENTAL  STUDY 


The  experimental  study  of  plasma  -was  conducted  in  quartz  chambers, 
circular  in  cross  section,  with  an  internal  diameter  of  a  =  0.6  cm. 
Cylindrical  chambers  40  cm  in  length  d  (Figure  2, a)  were  employed 
in  measurements  within  the  uniform  and  "corrugated"  magnetic  field 
(Sect.  3).  For  measurements  within  the  toroidal  magnetic  field 
(Sect.  4),  we  used  cylindrical  chambers  60  cm  in  length  d  bent  in 
the  form  of  a  circular  arc  with  a  radius  R  of  28  cm,  that  is, 
chambers  representing  part  of  a  torus  (Figure  2,b).  At  the  ends 
of  the  discharge  chambers  were  hot  tungsten  cathodes  and  anodes. 

A  magnetic  field  intensity  of  up  to  2600  Oe  was  created  using  a 
solenoid  consisting  of  individual  coils.  Magnetic  field  nonuni¬ 
formity  along  the  length  of  the  chamber  (along  its  axis)  did  not 
exceed  6%.  The  chambers  were  centered  relative  to  the  solenoid 
using  a  special  guide  tube  and  spacers. 

Prior  to  the  measurements,  the  experimental  chambers  were  filled 
with  spectrally  pure  helium  (pressure  p  =  0.02  -  1  mm  Hg).  The 
chambers  were  evacuated  over  a  period  of  40-50  hr  and  heated  at  a 
temperature  of  400-450°C  with  forced  heat.  The  chambers  were  then 
discharge-treated  and  filled.  The  pressure  of  the  gases  remaining 

in  the  system  detached  from  the  pumps  was  not  greater  than 

10-7  mm  Hg.  Filling  the  chambers  was  accomplished  in  a  number  of 
experiments  with  the  use  of  special  ampules  of  spectrally  pure 
gas  placed  in  the  discharge  chambers  prior  to  evacuation.  During 
the  measuring  procedure,  the  gas  pressure  was  monitored  with  the 
use  of  a  thermocouple  manometer  whose  graduation  was  preliminarily 
checked  using  a  McLeod  manometer. 

Plasma  was  created  within  the  chamber  by  means  of  high-voltage 
pulses  2(j.s  in  duration  and  at  a  repetition  frequency  of  30  Hz. 
Discharge  current  for  the  duration  of  the  pulse  amounted  to  0.5-1  A. 

The  concentration  of  charged  plasma  particles  was  then  measured 
during  the  period  between  discharge  pulses.  Measurements  were  per¬ 
formed  employing  the  resonator  method  (see  [3]>  for  example).  We 
placed  a  cylindrical  resonator  (internal  diameter  a^=2  cm,  length — 
1.8  cm)  with  apertures  in  the  ends  corresponding  to  the  size  of  the 
chamber  (Figure  2)  over  the  chamber  with  the  plasma.  The  resonator 
was  designed  to  generate  TMq^q  oscillations  within  the  3-centimeter 
wave  range.  In  the  case  of  this  type  of  oscillation,  the  electrical 
field  runs  parallel  with  the  axis  of  the  resonator.  The  loaded. Q- 
factor  for  the  resonator  with  the  quartz  chamber  was  1500-2000.  The 
connection  between  the  shift  in  resonator  resonance  frequency  &/, 
caused  by  the  introduction  of  the  plasma  into  the  resonator  ana 
the  average  concentration  of  electrons  with  respect  to  cross  sec¬ 
tion  n  is  given  by  the  relationship  (3] 

M _ L  A  JL.  n.  — 

/  ~  2  r  ’  Tp>  ~ 
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Figure  2.  1  -  discharge  chamber;  2  -  cathode;  3  -  anode;  4  - 

manometric  bulb;  5  -  measuring  resonator;  6  -  input  and  output 

waveguides;  7  -  solenoid  coil. 


where  /  -  field  frequency;  A ,  -  coefficient  of  the  form,  the  value 

of  which  is  calculated  with  the  following  approximated  relation¬ 
ship  with  the  selected  dimensions 


We  should  point  out  that  the  relationship  between  frequency  shift 
and  electron  concentration  does  not  depend  on  a  constant  magnetic 
field,  since  the  high-frequency  electrical  field  is  parallel  to 
the  magnetic  field.  The  effect  of  the  magnetic  field  on  the 
measurements  of  concentration  could  therefore  by  caused  only  by 
resonator  distortions.  The  effect  here  was  not  a  substantial  one, 
however,  since  the  measurements  were  performed  at  frequencies  far 
from  the  frequency  of  electronic  cyclotron  resonance  (in  contrast 
with  the  experiments  in  work  [!}). 

The  procedure  employed  in  measuring  the  frequency  shift  in  the 
case  of  the  resonator  with  the  plasma  was  analogous  to  that  em¬ 
ployed  in  work  [4],  and  we  will  not  devote  any  time  to  it.  A 


block  diagram  of  the  experimental  apparatus  employed  is  shown  in 
Figure  3.  The  accuracy  of  frequency  shift  measurements  was 
0.2  MHz.  The  minimum  measurable  concentration  was,  accordingly, 
in  the  neighborhood  of  2*10°  l/cm^. 

3.  RESULTS  OF  A  STUDY  OF  PLASMA  DECAY 
IN  A  UNIFORM  AND  CORRUGATED  MAGNETIC  FIELD 


Before  studying  the  effect  of  toroidal  nonuniformity  in  the  mag¬ 
netic  field  on  plasma  decay,  we  felt  that  we  should  first  under¬ 
take  a  study  of  plasma  decay  within  a  uniform  magnetic  field.  To 
conduct  these  experiments  we  employed  straight  cylindrical  cham¬ 
bers  (Figure  2, a)  filled  with  helium  at  a  pressure  of  0.04-0.4  mm  Hg. 


Figure  3.  GZ  -  activating  pulse 
generator;  GI  -  high-voltage 
pulse  generator;  BZ  -  delay  unit; 
GP  -  sawtooth  voltage  generator; 
KG  -  klystron  generator;  0  - 
oscillograph;  RB  -  discharge 
chamber;  R  -  resonator;  A  -  at¬ 
tenuator;  V  -  wave  meter;  D  - 
detector  section. 


Figure  4.  a )  ^  =  mh»  -  uniform 

magnetic  field  <»<«•/.>.  p  =  o.o4-om  Hg. 

1  -  data  from  the  present  study; 

2  -  data  from  work  [6],  p  =  0.35- 
-0.4  mm  Hg;  3  -  data  from  the 
present  study;  4  -  data  from 
work  [6];  b)  p  =  0.4  mm  Hg 


r~hHm 

«-»  = 
4 

~=m 

*.>*r 


a)- 

50%.  >)  p  i=  0.04-0.05 

.  .  .. 


In  the  case  of  charged  plasma  particle  concentrations*  lower  than 
5*10^-  I010cm"'5,  the  change  in  concentration  over  time  proved  to  be 


*We  refer  here  as  well  as  below  to  average  plasma  concentration 
with  respect  to  the  cross  section  of  the  chamber. 
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exponential.  This  is  evidence  of  the  fact  that  at  these  concen¬ 
trations  only  a  small  role  is  played  by  nonlinear  processes,  that 
is,  by  the  diffusion  connected  with  electron-ion  collisions  and 
the  recombination  (see  [5]).  The  time  constant  of  plasma  decay 
was  calculated  on  the  basis  of  the  slope  of  the  curves  describing 
the  dependence  of  the  logarithm  of  concentration  on  time 


1  _  rf(lrm)  (28) 

T  —  dt  '  V  ' 

Figure  4  presents  the  results  of  the  calculation  of  '  .  Plotted 

alone:  the  abscissa  is  the  value  —  (for  a  cylindrical  chamber 

4  Oj  * 

"#~T.8T  )*  which  in  the  case  of  the  diffusion  mechanism  of 

decay  should  be  equal  to  the  coefficient  of  lateral  diffusion. 

As  we  can  see  from  Figure  4,  the  results  from  our  calculation  of 
the  value  of  a  in  the  present  work  are  in  good  agreement  with 
the  results  presented  in  work  [6],  The  values  of  the  coefficient 

of  lateral  diffusion  prove  to  be  significantly  greater 

than  the  theoretical  values  (for  greater  detail,  see  [6]). 


We  also  tested  for  the  impact  on  plasma  decay  of  a  longitudinal 
"corrugated"  nonuniform  magnetic  field.  This  test  was  undertaken 
in  connection  with  the  fact  that  there  was  considerable  longitudi¬ 
nal  nonuniformity  (this  nonuniformity  reached  6$,  as  indicated 
above)  in  the  solenoid  comprising  individual  coils.  The  solenoid 
coils  were  combined  into  two  groups  in  order  to  create  a  controlled 
longitudinal  nonuniformity  in  the  magnetic  field.  The  current 
through  each  group  of  coils  was  varied  independently.  The  spatial 
period  of  the  nonuniform  field  (along  its  axis)  was  equal  to 
A 10  cm  .  During  the  measuring  process  the  cavity  resonator, 
which  was  used  to  determine  the  concentration  of  charged  plasma 
particles,  was  located  in  the  area  of  the  maximum  magnetic  field. 
The  results  of  measurements  of  the  decay  constant  within  a  range 

of  concentrations  of  2*10  -  5*10^  cm’^  are  presented  in  Figure  4. 

At  a  helium  pressure  of  * 0.4  mm  Hg,  the  time  of  longitudinal  ambi- 
polar  diffusion  over  a  length  equal  to  the  spatial  period  of  the 
nonuniformity  is  substantially  greater  than  the  time  of  lateral 
diffusion 


A* 


D 


15* 


(29) 


(in  estimation  we  assumed  On  =1400  cm^s  in  accordance  with  [7]j 
we  estimated  Chon  the  basis  of  the  experimental  data).  The  rate 
of  diffusion  had  therefore  to  be  calculated  by  the  magnetic  field 
at  the  point  where  the  measuring  resonator  was  located  (H=H  ). 
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The  relationship  between  — i  and  H  for  this  pressure  is  accord- 

t  max 

ingly  presented  in  Figure  4.  As  was  to  be  anticipated,  in  the 
case  of  a  different  degree  of  nonuniformity,  this  relationship 
proved  to  be  identical  to  the  one  established  in  the  case  of  the 
uniform  magnetic  field. 

If  longitudinal  diffusion  occurs  much  more  rapidly  than  lateral 
[if,  that  is,  the  inequality  the  inverse  of  (29)  is  observed], 
longitudinal  diffusion  in  the  region  of  the  minimum  magnetic  field 
with  subsequent  lateral  diffusion  is  the  shortest  "path  of  diffu¬ 
sion"  from  out  of  the  region  of  the  maximum  magnetic  field.  Eif- 
fusion  time  should  accordingly  be  near  the  time  of  lateral  diffu¬ 
sion  in  the  case  of  a  magnetic  field  equal  to  the  minimum.  The 
conditions  of  the  experiments  performed  at  a  helium  pressure  of 
0.04-0.05  mm  Hg  approach  this  case.  As  we  can  see  from  Figure  4, 
the  relationship  between  and  Hmin  in  the  case  of  the  nonuniform 

magnetic  field  at  this  pressure  is  found  to  be  the  same  as  that 
occurring  in  the  case  of  the  uniform  field. 

The  results  presented  indicate  that  plasma  decay  time  in  a  magnetic 
field  with  longitudinal  nonuniformity  is  found  to  be  the  interval 
between  the  times  of  lateral  diffusion  with  the  maximum  and  minimum 
magnetic  field. 

4.  RESULTS  OF  A  STUDY  OF  PLASMA  DECAY  IN  A  TOROIDAL  MAGNET TO  FIELD 

The  discharge  chambers  shown  in  Figure  2,b  were  used  to  study  plasma 
decay  within  a  toroidal  magnetic  field.  The  chambers  were  filled 
with  helium  at  a  pressure  of  0.02-0.9  mm  Hg. 

Measurements  indicated  that  under  certain  conditions  (at  pressures 
of  0.05-0. 5  mm  Hg  and  with  magnetic  fields  greater  than  1000-2000  Oe) 
the  change  over  time  in  charged  particle  concentration  differs  from 

the  exponential  within  the  range  of  values  2*10  -  10  cm*,  in 
which  nonlinear  processes  are  unimportant,  A  distinct  disconti¬ 
nuity  was  observed  Ln_this  regard  on  the  curves  describing  the 
relationship  betweeirm  and  t.  Characteristic  curves  of  this  type 
are  shown  in  Figure  5.  In  all  cases  in  which  a  discontinuity  was 
observed,  the  value  of  the  decay  constant,  as  calculated  with  re¬ 
spect  to  the  slope  of  the  curves  in  the  initial  region  (at  concen¬ 
trations  of  2*10^-  1010  1/crn-^),  was  close  to  that  of  the  decay 
constant  within  the  uniform  magnetic  field  (in  straight  chambers). 
Results  from  the  calculation  of  the  plasma  decay  constant  with 
respect  to  the  slope  of  the  curves  In n=/(t)  at  the  termination  of 

decay  (with  concentrations  lower  than  10^  cm  ^)  are  presented  in 

Figure  6,  The  relationship  between  —  and  H  has  a  characteristic 
minimum.  It  was  possible  to  observe  this  minimum  with  magnetic 
fields  weaker  than  2600  Oe  for  helium  pressures  lower  than  0.15  mm  Hg 
At  a  helium  pressure  of  0.4  mm  Hg  we  observe  only  a  limited  slowing 

down  in  the  drop  of  curves  -J-  ( H )  as  compared  with  the  curves 


obtained  for  the  uniform  magnetic  field  (in  the  straight  chamber). 
At  a  pressure  of  0.9  mm  Hg,  the  relationship  between  t  and  H  is 
found  to  be  virtually  the  same  as  that  established  for  the  uni¬ 
form  field.  ' 


Let  us  now  compare  the  experimental  data  with  the  results  obtained 
from  theoretical  analysis.  Let  us  note  first  of  all  that  under  the 
conditions  prevailing  in  the  experiments  which  were  performed,  lon¬ 
gitudinal  diffusion  occurs  much  more  slowly  than  lateral  diffusion 


The  characteristics  of  plasma  decay  in  the  chamber 


representing  part  of  a  torus  (Figure  2,b)  should  accordingly  be 
identical  to  those  in  a  closed  torus. 


The  shape  of  the  curves  of  plasma  decay  similar  to  those  shown  in 
Figure  5  may  be  easily  explained  if  we  assume  that  at  the  beginning 
of  the  decay  process,  the  distribution  of  concentration  over  the 
cross  section  of  the  chamber  is  symmetrical  and  close  to  the  diffu¬ 
sion  [concentration — TR],  This  would  appear  to  be  a  plausible  as¬ 
sumption,  since  the  nonlinear  processes  of  removing  the  charged 
particles,  which  at  the  onset  of  decay  may  be  more  effective  than 
toroidal  drift,  lead  to  a  " symmetrization"  of  the  distribution  of 
concentrations.  In  the  case  of  a  symmetrical  distribution,  the 
rate  of  plasma  decay  should,  according  to  the  theory,  be  identical 
with  that  in  the  uniform  magnetic  field  (p,  1047).  This  is  pre¬ 
cisely  what  is  to  be  observed  in  the  experiment.  Only  after  a 
time  equal  to  several  "drift  times"  should  we  be  able  to  establish 
the  nonsymmetrical  distribution  of  concentrations  characteristic 
of  the  toroidal  chamber  (24).  1 


i 
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Let  us  now  compare  relationship  (25),  which  defines  the  constant  of 
plasma  decay  after  the  distribution  of  concentrations  has  been  es¬ 
tablished,  with  the  experimental  results  shown  in  Figure  5.  We 
should  keep  in  mind  that  in  the  experiment,  the  value  of  the  co¬ 
efficient  of  lateral  diffusion  in  the  uniform  magnetic  field  is 


Figure  7,  Experimental  data:  1  -  p  =  0,12;  2  -  0,055;  3  -  0,04; 
4  -  0.025  Hg;  Solid  line  represents  theoretical  curve. 
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found  to  be  substantially  greater  than  that  arrived  at  theoretically 
(see  p.  1051  and  [6]).  In  making  our  comparison,  therefore,  we  will 
substitute  in  (25)  the  experimentally  calculated  values  of  the  co¬ 
efficient  of  lateral  diffusion,* 

In  the  accompanying  table  we  compare  the  experimental  and  theoret¬ 
ical  values  of  the  maximum  decay  constant  and  the  corresponding 
values  of  the  coefficient  of  lateral  diffusion.  The  theoretical 
values  xm  and  D0  are  calculated  in  accordance  with  the  formula 
in  (25).  In  accordance  with  the  data  presented  in  [7]#  we  have 
assumed  the  value  of  thf*  coefficient  of  longitudinal  diffusion 
here  to  be  equal  to  f^it=540/p  cm2/s  (p  in  mm  Hg).  Experimental 


1 , 

mm  pt.  CT. 

**,,  mu—m.Xf  | 

TIAtkHI 

3 

CIM4 

1  •BCO«pHM«a~ 

-  TIAkaWt 
\  3  IM44UI 

mp«T8W 

\t) 

0.025 

100 

80 

70 

95 

0.04 

130 

130 

70 

60 

0.055 

160 

180 

35 

43 

0.12 

280 

380 

20 

20 

0.4 

>400 

1300 

<30 

6 

Key  to  Table:  1  -  mm  Hg;  2  -  ps;  5  -  experimental  values;  4  -  the- 

p 

oretical  values;  5  -  cm  /s. 


values  of  T.  were  calculated  using  the  curves  in  Figure  6.  The 
value  of  the  magnetic  field  with  which  the  decay  constant  is  maxi¬ 
mum  was  calculated  on  the  basis  of  these  curves  as  well.  The  con¬ 
nection  between  the  magnetic  field  and  the  coefficient  of  lateral 
diffusion  was  established  on  the  basis  of  the  results  of  measure¬ 
ments  performed  in  the  uniform  magnetic  field. 

We  can  see  from  the  table  that  at  a  helium  pressure  of  0.02-0.12 
mm  Hg,  the  experimental  and  theoretical  values  D0  are  in  good 
agreement  with  one  another.  The  case  in  which  at  pressures  of 
0.4  and  0. 9  mm  Hg  a  minimum  on  curves  JL(//)  was  not  observed  (with 

H  <2600  Oe)  is  also  in  agreement  with  the  theory. 


Let  us  now  compare  the  path  of  curves  ~  ( H )  for  various  pressures. 
Figure  7  shows  the  relationship  between  J_and  pk  in  relative 


*It  is  assumed  in  the  case  of  this  implicit  example  that  the  theo¬ 
retical  expression  for  the  coefficient  of  lateral  diffusion  is  cor¬ 
rect,  while  deviations  of  the  experimental  values  from  the  theoret¬ 
ical  are  connected  with  an  abnormally  high  frequency  of  collisions. 
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unit .  This  relationship  is  constructed  by  employing  the  curves 
in  Figure  6  and  the  experimental  curves  describing  the  relation¬ 
ship  between  D.  and  H  within  a  uniform  magnetic  field.  The  ex¬ 
perimental  data  obtained  at  various  helium  pressures  all  fit  the 
same  curve  in  accordance  with  (25).  The  values  of  *  on  the  curve 
differ  from  the  theoretical  values  by  no  more  than  a  factor  of  1*5* 

The  experiments  which  have  been  performed  thus  indicate  that  the 
theory  set  forth  above  correctly  describes  the  impact  of  toroidal 
drift  on  plasma  decay. 

The  authors  express  their  profound  gratitude  to  B.  B.  Kadomtsev 
and  V.  I.  Perel 1  for  helpful  discussions  of  the  results  of  their 
work. 
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